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1. Introduction 

The exclusion process is one of the most basic stochastic interacting particle 
processes. It consists of particles executing independent random walks on a graph, 
subject to the exclusion interaction that prevents two particles from occupying the 
same vertex. In this paper we look at the version where the underlying graph is 
the integer lattice Z with nearest-neighbor bonds. The particles take only nearest- 
neighbor steps to the right, and a jump is permitted only if the next site to the 
right is vacant. In the spatially homogeneous case jumps across all nearest-neighbor 
bonds happen at exponential rate one. This model goes by the name TASEP, or 
totally asymmetric simple exclusion process. 

In the early 1990's Janowsky and Lebowitz [JL1, JL2] introduced an interesting 
variant of TASEP. Jumps across a specific bond, say {0, 1}, happen at a rate r 
strictly lower than the rate 1 everywhere else. Imagine for example a tollbooth on 
a single lane highway that slows down the flow of cars. For the slow bond model 
one investigates the same questions as for homogeneous TASEP, such as large scale 
behavior and invariant distributions. Simulations and partial results are available, 
but crucial results remain open. 

In the present paper we prove a hydrodynamic limit for the slow bond model. In 
one sense our result is complete and general. It is valid for all values r G (0, 1) of the 
slow rate. We need no extra assumptions on the initial distribution of the particles, 
only that the law of large numbers for the density profile be valid at time zero. 
But the hydrodynamic limit cannot make the most interesting distinction, namely 
whether the slow bond disturbs hydrodynamic profiles for all values r < 1. This 
question can be resolved only through sharper control of particle-level evolution. 

The macroscopic effect of the slow bond can be described like this: Corresponding 
to the microscopic rate profile (rate identically 1 at all sites except 0, and rate 
equal to r at 0) there is a macroscopic rate profile X(x) defined for x G R such 
that X(x) = 1 for x ^ 0, and A(0) G (r,l\. The macroscopic evolution of the 
particle density is the solution of an optimal control problem whose running cost 
depends on the space variable x through the function X(x). If A(0) = 1, there is no 
visible macroscopic disturbance from the slow bond. The interesting open question 
becomes: Is A(0) < 1 for all r < 1? 

We give a precise definition of the quantity A(0) in terms of a last-passage growth 
model whose association with TASEP is well-known. In this formulation the open 
problem can be cast in a percolation-type framework, and perhaps made amenable 
to other techniques, such as the random matrix methods that have recently been 
applied to growth models with great success. 

To prove the result we use the variational coupling method initiated in [Sel] 
which, when it applies, gives laws of large numbers without any knowledge of in- 
variant distributions. This approach allows us to work directly on the slow bond 
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model, without intermediate approximations with more tractable models. The 
proof extends to the situation with multiple slow bonds, and also works for more 
general exclusion processes that admit K particles per site. 

In addition to the Janowsky-Lebowitz papers mentioned above, there is a hand- 
ful of other related work. Covert and Rezakhanlou [CR] derived a bound for the 
critical value of the slow rate by approximating the slow bond model with an exclu- 
sion process whose rates vary more regularly in space. Liggett's [Li] new monograph 
discusses the slow bond model and proves bounds for the critical value of the rate. 
As for Janowsky and Lebowitz, the approach is through finite systems with open 
boundaries, with system size tending to infinity. The case of a zero-range process 
(ZRP) with a slow site was treated by Landim [La]. In ZRP the particles accumulate 
at the slow site and produce a point mass in the hydrodynamic profile. The ZRP 
is a more tractable model because it retains product-form invariant distributions 
even when rates lose spatial homogeneity. 

Our paper is organized as follows. The main result is the hydrodynamic limit 
Thm. 2.2. On the way to it we define the macroscopic rate A(0) in terms of the 
growth model, and state some bounds for it in Thm. 2.1. As corollaries to Thm. 
2.2 we compute the macroscopic profiles that evolve from constant initial profiles, 
characterize the macroscopically invariant profiles in the range [p*, 1 — p*] between 
the lower and upper critical densities, and prove one property of invariant measures. 
These results are in Sect. 2. The remainder of the paper contains the proofs. 

2. Results 

The process operates according to these rules: Indistinguishable particles occupy 
the sites of the one-dimensional integer lattice Z. Each site has at most one particle, 
so the state of the process is described by the occupation numbers rj = (i]i)i e z : 
where r\i = 1 if site i is occupied, and if site i is empty. Particles take nearest- 
neighbor steps to the right, subject to the exclusion rule, at exponential rate 1, with 
one exception: jumps from site to site 1 happen at rate r. The rate r is a fixed 
constant in the range (0, 1]. On the compact state space {0, 1} Z , the dynamics has 
the infinitesimal generator 

(2.1) Lf( v ) = J2v l (i-m + i)[f(v hl+1 )-f(v)] + r V0 (l - mMv^ 1 ) - f( V )} , 

where rj l,l+1 = r\ — 5i + 5i + i denotes the configuration that results after a single 
particle jumps from i to i + 1. When r = 1 this is the generator of TASEP, the 
totally asymmetric simple exclusion process. 

We prove a hydrodynamic limit for this process in the usual Euler scale, for all 
values of r. As is well-known, when r = 1 the macroscopic particle density p(x,t) 
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obeys the scalar conservation law 

(2.2) p t + f (p) x = 

with current 

(2-3) / (p) = p(l-p). 

Starting with Rost [Ro] in 1981, this hydrodynamic limit of the space- homogeneous 
TASEP has gone through many stages of generalization and refinement. See Ch. 8 
of [KL], Part III of [Li] and their notes and references. 

The slow bond restricts the range of admissible currents. In TASEP the range 
of currents is [0,1/4], and the maximal current 1/4 = /q(1/2) occurs at density 
p = 1/2. In the slow bond model the maximal current is some value J max ( r ) < 1/4. 
We know rigorously that J max .(r) < 1/4 if r is small enough. When this happens, 
the densities p around 1/2 for which fo(p) > J m ax(r) become inadmissible. This 
we can prove on the hydrodynamic scale. 

Let n(r) = J m ax(^) _1 be the reciprocal of the maximal rate. Thus k(t) is a 
nonincreasing function of r G (0, 1], with values in the range [4, oo) and k(1) = 4. 
Since we do not know the exact dependence of n(r) on r, we shall give it a precise 
definition through a well-known growth model. 

2.1. Definition of n(r) 

Consider the following last-passage growth model on the first quadrant of the 
plane. For e N 2 , let Y iy j be i.i.d. exponentially distributed random variables 
with common expectation E[Yi j] = 1. Let the deterministic weights Wi j be given 
by 

(2-4) w hJ = < . . 

{ 1/r, i = j- 

For n = 1, 2, 3, . . . , let the passage time of site (n, n) be 

2n-l 

(2.5) T n>n = max \ w i h ,j h Yik,ik > 

7T ' * 



k=l 



where the maximum is over paths n = {(1, 1) = (h, ji), («2, J2), • • • , (*2n-i> hn-i) = 
(n, n)} that take steps only to the right and up: for each k, 



(2.6) 



(h,jk) - (ik-i,3k-i) = (1,0) or (0, 1). 
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An obvious superadditivity holds for the random variables {T n n }, and one can 
derive moment bounds (see for example Thm. 6.3 in [GW], or Prop. 5.1 in [Se3]) 
sufficient for Kingman's [Ki] subadditive ergodic theorem, to obtain a strong law 
of large numbers: There exists a constant «(r) such that 

1 

(2.7) lim -T nn = n{r) a.s. 

n^oo fl 

This limit is taken as the definition of n(r). It is known that k(1) = 4. This is a 
special case of the full interface result: for (x, y) G R+ 

1 2 

(2.8) lim -T [nx] [ny] = U/x + y/y) a.s. in the case r = 1. 

n^oo n 

Proofs of (2.8) can be found in [Ro] and [Sel]. 
Currently we have this information on n(r): 

Theorem 2.1. The function n{r), < r < 1, is continuous and nonincreasing. It 
satisfies these bounds: 

< kOti) — k{t2) < — — for ri < r 2 , 

and 

(2.9) max(4, - + - + 2(1 + r) \ < ft ( r ) < 3 + - . 
y 1 \ ' 2 2r(l +r) J " V ; " r 

The proof shows that neither bound in (2.9) is optimal. Janowsky and Lebowitz 
[JL2] give the bound J max (^) < ^(1 — f) for r < 1/2. (For the proof, see p. 277 in 
[Li].) For k{t) this gives 

(2.10) K{r) -^hr) fOT ^V2. 

This improves (2.9) in the range r G [(vT7 — l)/8, 1/2]. However, the definition of 
J max (r) in [JL2] and [Li] is not the same as ours, which is tied to the hydrodynamic 
limit. In [JL2, Li], J max (0 is the limiting stationary current of a finite system 
with open boundaries, as the size of the system tends to infinity. Of course the 
two quantities ought to be the same. Once optimal bounds are found with some 
approach, it will be of interest to verify that the different definitions are the same. 
Let us define 

r* = inf{r G (0, 1] : «(r) = 4} . 
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The interesting open problem is whether r* = 1. In other words, does the macro- 
scopic passage time n(r) rise strictly above «(1) = 4 as soon as r < 1? Equivalently, 
is there a forbidden range of densities around 1/2 as soon as r < 1? Simulations 
by Janowsky and Lebowitz [JL2] suggest that such is the case. The lower bound 
in (2.9) gives r* > (-v/41 - 3)/8 « 0.425, while (2.10) gives r* > 1/2. 

2.2. The hydrodynamic limit 

Let po be a given measurable function on R such that < po(x) < 1 for all 
x G R. This is the initial macroscopic profile. Let vq be the antiderivative of po 
defined by 

(2.11) uo(0) = 0, v (b)-v (a)= I p (x)dx . 

J a 

Define the "macroscopic rate profile" X(x) as 

r 1, Xy^O 

(2.12) X(x) = I 4 
where n(r) is defined by (2.7). For x G R, let 



x = 0, 



g (x) = sup {/o(p) ~xp} 

0< P <1 

denote the Legendre conjugate of fo- It is given by 

-x , x < — 1 

(2.13) 9o(x) = { (l/4)(l-x) 2 , -1<x<1 

0, x > 1 . 

For x G R and t > 0, define v(x, 0) = fo(a^), and for t > 0, 

(2.14) = supiuoMO)) - / A(w(a))ff f w ( ^ V|- 

tu(.) I Jo \X{w{s))J J 

The supremum is over piecewise C 1 paths w : [0,£] — > R that satisfy w(t) = x. In 
Section 5.1 we give a formula for the path w(-) that minimizes the integral part 
inside the braces for a given initial point w(0) = q. 

For each fixed t, v(- , t) is a Lipschitz function. Its cc-derivative 

(2.15) p(x,t) = -^-v(x,t) 



7 



represents the macroscopic density profile of the particles. 

If X(x) = 1 for all x, (2.14)-(2.15) give the entropy solution of (2.2) with initial 
data pq. This is the density profile of space-homogeneous TASEP. Formula (2.14) 
expresses the sense in which X(x) can be regarded as the macroscopic rate profile. 
From (2.9) we get 

4r 

A <°> * FTT 

This implies that A(0) > r for r < 1. In other words, the microscopic averaging 
does not simply reproduce the rate r at the macroscopic level. 

For the hydrodynamic limit, assume that we have constructed a sequence of 
exclusion processes i] n (t) = (77™ (t) : i G Z) for times t > 0, where n = 1, 2, 3, . . . is 
the index of the sequence. The initial distributions of the processes are arbitrary, 
subject to the condition that this weak law of large numbers is valid: 
for all a < b and e > 0, 

[nb] 

„ . > =0. 

i=[na] + l 

We wrote P n for the probability measure on the probability space of the process 
rf 1 . Let J™(t) denote the number of particles that have made the jump from site i 
to i + 1 in the time interval [0, £], in the process r? n (-). 

Theorem 2.2. Let t/ie slow rate r be any number in (0,1]. Under assumption 

(2.16) , these weak laws of large numbers hold at macroscopic times t > 0: for all 
real numbers a < b and e > 0, 

(2.17) lim P n {\n- 1 Jf na] (nt)-(v (a)-v(a,t))\>e} = 0, 



(2.16) 



r i [nb] r b i 

lim P n { - V <(0)- / ,9 (2;)cte 



and 



,- | i nh \ rb 

lim P n { - V r]?(nt)- p{x,t)dx 
n ^°° { n ,_f—' J a 



> e I = 



i=[na]+l 

where v(x,t) is defined by (2.14) and p(x,t) = v x {x,t). 

Let us derive some corollaries of Theorem 2.2. When A(0) < 1 there is a critical 
density 

(2-18) p* = \ - ^ v / l 3 A(bT, 

determined by the current condition fo(p*) = A(0)/4. The slow bond disturbs the 
hydrodynamic profile only for densities in the range p G (p*, 1— p*). When A(0) = 1 
[or, equivalently, the passage time k(t) =4], p* = 1/2 and the interval (p*, 1 — p*) 
is empty. 

As a first application of Thm. 2.2 we derive the macroscopic evolution of constant 
initial profiles. 
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Corollary 2.1. Assume the initial distribution of the particles is chosen so that 
(2.16) is satisfied for the constant profile po(x) = p. Then the macroscopic limits 
in Thm. 2.2 are described as follows: 

Case 1: pG [0, p*] U [1 — p*, 1] . Then v(x,t) = px — tp(l — p) and macros copically 
we see a constant density p(x,t) = p, exactly as for the homogeneous TASEP. 

Case 2: A(0) < 1 and p G (p*, 1 — p*). The behavior changes around the blockage 
at the origin: 



Case 1 is still valid for x outside the range in (2.19). Correspondingly, the density 
profile has constant segments of upper and lower critical densities around the origin: 



The reason for this behavior is that the maximal current permitted by the blocked 
system is A(0)/4 = p*(l — p*), while the unblocked TASEP has current fo(p) = 
p(l — p). The blockage does not disturb the system unless the system tries to 
transport particles at a current above A(0)/4. This happens if p(l — p) > A(0)/4, 
which is equivalent to p G (p*, 1 — p*). 

Corollary 2.2. Suppose A(0) < 1 so that p* G (0,1/2). Let po be a macroscopic 
profile that is piecewise continuous in each bounded interval, and satisfies p* < 
Po(x) < 1 — p* . Suppose po is invariant under the macroscopic dynamics, so that 
p(x,t) = po(x) for a.e. x, for all times t. Then Lebesgue- almost everywhere po(x) 
G {p*, 1 — p*}, and po takes 0-3 jumps according to these restrictions: a jump from 
1 — p* to p* can occur only at the origin, while jumps from p* to 1 — p* can occur 
at any x G R. 

We can control po only almost everywhere because a macroscopic density profile 
p(- ,t) is determined only through its integral v (• ,t). The upward jumps from p* 
to 1 — p* are the usual entropy shocks of equation (2.2). The slow bond permits a 
downward jump from 1 — p* to p* at the origin, which violates the entropy condition 
for the concave current fo- 

By Cor. 2.2 there cannot be a macroscopically invariant profile strictly in the 
range (p*, 1 — p*). Consequently, the process cannot have an invariant probability 
distribution for which a macroscopic profile exists in the sense of (2.16) and lies in 
the range (p*, 1 — p*). The interesting open problem is to prove the existence of 
an invariant probability measure that corresponds to the macroscopically invariant 
non-entropy shock profile po(x) = (1 — p*)l{x < 0} + p*l{x > 0}. About invariant 
measures we have this to say: 



(2.19) 




(2.20) 




9 



Corollary 2.3. For any valuer of the slow rate, and for any p E [0, p*]U[l— p*, 1], 
there exists an invariant distribution p for the process such that 



(2.21) ^ = 1,^+1 = 0} = 



r 1 p(l — p) , i = 
p(l-p), i + 0. 



In particular, suppose the initial distribution of the process is the Bernoulli dis- 
tribution «o with marginals 

cto{Vi = 1} = 1 ~~ P* f or i < 0; an d a o{Vi = 1} = P* f or i > 0- 

Let at denote the distribution at time t > 0. Then any limit point p of the time 
averages t~ l a s ds satisfies (2.21) with p = p* . 

2.3. Remarks and extensions 

K -exclusion and multiple slow bonds. Our proof uses the method of [Se3]. As 
in that paper, we could prove Thm. 2.2 also for generalized exclusion processes 
that allow K particles per site, instead of just one particle. The result would be 
qualitatively the same, but with a different growth model in Sect. 2.1 and corre- 
spondingly different A(0). With suitable large deviation estimates, as in [Se3], we 
envision that the weak law of large numbers in Thm. 2.2. could be strengthened 
to a strong law of large numbers (convergence for almost every realization of the 
microscopic evolution) . 

A natural extension is to consider finitely many macroscopically separated slow 
bonds. Suppose that for the nth process n n , jumps from site [ occur at rate 
rfc for some set x\ < • • • < x m of macroscopic space points. Thm. 2.2 remains 
unchanged. The change appears in the definition (2.12) of the macroscopic rate 
profile: the function A(-) would take on the values X(xk) = A/n{rk), and X(x) = 1 
elsewhere. 

Hamilton- Jacobi equations. (2.14) describes v(x,t) as the value function of an 
optimal control problem with running cost X(w(s))go(w' (s) / X(w(s))) . If the func- 
tion A had sufficient regularity, standard theory would imply that v(x,i) is the 
unique viscosity solution of the Cauchy problem 

(2.22) v t + X(x)f (v x ) = , v(- , 0) = ^- 

(See Ch. 10 in [Ev].) One can directly check that the solutions in Cases 1 and 2 
of Cor. 2.1 satisfy (2.22) except at jump points. But there is currently no theory 
about existence and uniqueness of solutions of Hamilton- Jacobi equations with a 
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discontinuity of the type that X(x) has at 0. The presently available approaches to 
discontinuous equations involve smoothing out the discontinuity with a mollifier, 
see [Os] and its references. But once A has been convolved with a mollifier, the 
discontinuity at is lost. So a different approach is required for our problem. 

The growth model. The connection between TASEP and the homogeneous ver- 
sion of the growth model of Section 2.1 goes back to Rost's 1981 paper [Ro]. The 
last-passage formulation was first utilized for hydrodynamic limits and large devi- 
ations in [Sel, Se2]. The connection between exclusion and the growth model will 
appear in Section 4. 

An additional motivation for the definition of «(r) in terms of the path model 
is this: analysis of models of this type with tools from combinatorics and random 
matrix theory has recently led to elegant exact calculations of limits and fluctua- 
tions, in the work of Baik, Deift, Johansson, and Rains. The homogeneous version 
of this particular growth model is treated in [Jo] . 

The hydrodynamic result of Thm. 2.2 implies a shape result for the last-passage 
growth model. We will not write down the details of such a conversion, but only 
note one point: Suppose the diagonal defect does not pass through the origin, but 
instead through the macroscopic point (0, u) for some u > 0. To capture this, 
redefine the weights in (2.4) by w iji+ [ nu ] = r _1 , and for other points Wij = 1. 
Now the weights change with n, to keep the defect line at the macroscopic line 
y = x + u. The limiting shape for this model follows (2.8) part of the way, but 
develops a kink around the point where the defect line passes through the interface, 
and the interface is no longer convex. This conclusion can be worked out from the 
explicit formulas of Section 5.1, via the mapping in Lemma 4.2. 

Relation with the Covert-Rezakhanlou result. The definition of «(r) through 
(2.4)-(2.7) suggests an immediate bound. If we increase all weights in (2.4) to 
Wij = 1/r and use Rost's result (2.8), we get the upper bound K,(r) < 4/r. By 
(2.12) this implies A(0) > r, and then from (2.18), p* > (1/2) - (l/2)y/T=r. 
Equivalently, the blockage does not disturb a profile at density p if r > 4p(l — p). 
This is the bound of Covert and Rezakhanlou [CR] . 

To obtain their bound [CR] approximated the blocked TASEP with a system 
whose jump rates depend on spatial location through a continuous function X(x), 
so that in the nth process jumps from site i occur at rate X(i/n). 

The Rost picture shows why a continuous, one-sided macroscopic approximation 
cannot get any closer than the upper bound n(r) < 4/r. Let X(x, y) be a continuous 
function that satisfies X(x,y) < 1 and \(x,x) = r. In (2.5) take new weights 
Wij = X(i/n, j/n) -1 , and compute k(r) as the limit in (2.7) with the weights Wij. 
Since Wij < Wij, we have k(t) < R(r). Given any e > 0, choose 6 > so that 
X(x,y) < r + e in a strip of width 5 around the diagonal A = {(x,x) : < x < 
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1}. In the homogeneous problem (weights = 1), the strict concavity of the limit 
(2.8) implies that, macroscopically, the diagonal is the unique maximizing path for 
(x,y) = (1,1). Thus with high probability, the n5-strip around the microscopic 
diagonal : 1 < i < n} contains a path n such that Yij > 4n + o(n). The 

weights in the n5-strip satisfy Wij > (r + e) -1 , so it follows that R{r) > 4(r + e)~ 1 . 
Since e > was arbitrary, k(r) > 4/r, and we see that the upper bound 4/r is the 
best we can get by continuously bounding the rates from below. 

3. Proof of Theorem 2.1 

Let us write T^ n for the passage time in (2.5) to indicate dependence on r. Let 

< T\ < T2 < 1- Let 7T be the path that gives the maximum in (2.5) for T^ n: so 
that 

^ = £y u+ (1-i) £ Y u . 

u£7T ^ ' u67rnA 

We wrote u = for an integer site on the plane and A = {(x, x) : x G R} is the 
diagonal. Since we are maximizing passage times of paths, 

T^>$>+(^-l) £ Yu, 

u£7r ^ ' u67rnA 

and we get 

T n,n ~ T n% < ( ~ ~ — J £ Y i,i • 

Divide by n and let n — > oo to get the first inequality in Thm. 2.1. 

Set T2 = 1 and use k(1) = 4 to get the upper bound in (2.9). The value n(l) = 4 
is part of the lower bound because k(-) is a nondecreasing function. 

To get the other part of the lower bound, consider this particular path ir n from 
(1, 1) to (n, n). For each % = 1, . . . , n — 1 choose the larger one of Y i+ i }i and Y iji+ i. 
Since £7[Y' V7"] =3/2 for two independent rate one exponentials, this contributes 
(approximately) 3n/2 to the sum. For each pair {z,z + 1}, the above choice led to 
one of two situations: 

Case 1. Two points on the same side of the diagonal: (i, i + 1) and (i + 1, i + 2), 
or (z + 1, z) and (i + 2, z + 1). 

Case Two points on opposite sides of the diagonal: + and (i + 2, z + 1), 
or (z + 1, z) and (z + 1, z + 2). 

To complete the path it we pick one site between each pair chosen above. In Case 

1 we can choose the larger one of the diagonal r -1 Y i+ i 5 j + i and an off-diagonal value, 
Y^i_i_2 or Yi + 2,i depending on the subcase of Case 1. The expected contribution 
is E[Y' V (r -1 !"")] = (r 2 + r + l)r _1 (l + r) -1 . In Case 2 there is no choice: 
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the path must go through (i + l,i + 1) and take the diagonal value r~ 1 Y i+ i^ i+ i. 
Cases 1 and 2 are equally likely, so this second step contributes approximately 
(n/2) ■ (r 2 + r + l)r _1 (l + r) -1 + (n/2) ■ r -1 . Adding the contributions from the 
two steps gives the lower bound in (2.9). 

4. Proof of Theorem 2.2 

4-1. Construction of the process and the variational coupling 

We follow the construction in Section 4 of [Se3], and only outline it here. We 
construct a process z(t) = (zi(t) : i E Z) of labeled particles that move on Z. The 
location of the ith particle at time t is Zi(t), and these satisfy 

(4.1) 0<z i+1 (t)-Zi(t)<l. 

For the dynamics, let {T>i} be a collection of mutually independent Poisson jump 
time processes on the time line (0, oo). T>q has rate r, and all other X>j have rate 1. 
In the graphical construction, Zi attempts a jump one step to the left at epochs of 
T>i. The jump is executed if it does not violate (4.1). This happens independently 
for all i. 

Once the z(t) process is constructed, the exclusion process 77(f) is defined by 

(4.2) m(t) = Zi (t) - Zi-^t) . 

The ^-particles keep track of the current of the exclusion process rj(-). The number 
of r/-particles that have left site i in time [0, t] is given by 

(4.3) Ji(t) = Zi (0) - Zi (t). 

Assume now that the process z(-) has been constructed on some probability space 
that supports the initial configuration z(0) = (^(0)), and the Poisson processes 
{T>i} that are independent of (^(0)). We define a family {w k : k E Z} of auxiliary 
processes on this same probability space. Each w k (t) = (wf(t) : i E Z) is an 
exclusion process just like z{t), so (4.1) is in force for all i. Initially 



(4.4) w^O) 



z k {0), z>0 
z k (0) + i, i<0 



and dynamically 

(4.5) w\ attempts to jump to w\ — 1 at the epochs of V i+k . 

The usefulness of the family of processes {w k } lies in this fact: 
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Lemma 4.1. For all i G Z arid t > 0, 

(4.6) = supu;?_ fc (t) a.s. 



This lemma is proved as Lemma 4.2 in [Se3] so we will not repeat the proof here. 
This is the "variational coupling" that is the key to our proof. 

Thinking of the process w k {t) as giving the height of an interface over the sites i, 
we want to normalize it to start at height zero, and also to advance in the increasing 
coordinate direction. Hence define a new family of processes {£ fc } by 

(4.7) tf(t) = z k (0)-w>°(t) forzeZ, t > 0. 
Now we can write (4.6) as 

(4.8) zi(t) = sup{z*(0) - £{_*(*)} . 

kez 

The virtue of (4.8) is that the effect of the initial condition (zfc(0)) has been sep- 
arated from the effect of the Poisson jump times {T>i}. The process £ fc does not 
depend on z k (0), and depends on the superscript k only through a translation of 
the indexing of {T>i}. Initially 



(4.9) tf(0) 



0, i > 
-i, i < 



The dynamical rule for the process is that 

(4.10) £f jumps to + 1 at epochs of T> i+k , 



provided the inequalities 

(4.H) &<&-i and ef<ef+i + 1 



are not violated. Notice that £_ fc jumps at rate r, while other £f jump at rate 1. 

We can now outline the strategy for proving Thm. 2.2. Given the initial con- 
figurations r] n (0) = (r}f(0) : i G Z) that appear in hypothesis (2.16), define initial 
configurations z n (0) = (zf (0) : i G Z) so that Zq(0) = an< ^ (4-2) holds at time 
i = 0. Then hypothesis (2.16) implies that 

(4.12) lim n~ 1 z?! lol (0) = v (q) in probability 



for all g G R, with v defined by (2.11). 



14 



Construct the processes z n (t) as indicated above, and define the exclusion pro- 
cesses rj n (t) by (4.2). Define v(x,t) by (2.14). From (4.2)-(4.3) we see that both 
limits of Thm. 2.2 follow from proving that for all x G R and t > 0, 

(4.13) lim n~ 1 ZL nx] (nt) = v{x, t) in probability. 

n — >oo L ' 

Now rewrite (4.8) with the correct scaling: 

(4-14) n-'z^nt) = sup {n" 1 ^, (0) - n' 1 ^^(nt)) . 

In this formula each process z n (-) is defined on a probability space that supports the 
initial configuration z n (0) and the Poisson processes {T>i}. On each such probability 
space we define the processes {£ fc (-)} as functions of according to (4.9)-(4.11). 
The proof of (4.13) is now to show that the right-hand side of (4.14) converges to 
the right-hand side of (2.14). The first step is to prove a limit for the £-term. 

4-2. Limit for £ 

First some definitions. The meaning of these notions will be explained below. 
Let 

V = {(x,y)GR 2 :y>0,x>-y}. 

For (x, y) G V, let 

(4.15) Ko(x,y) = + Vv) 2 - 

Let x(s) = (xi(s),X2(s)) denote a path in R 2 , defined on some interval of s-values. 
For (x, y) G V and q G R let 

(4.16) r«(x, y) = sup( f J 0{ \'[ S) \ ds : x(-) G X(x, y) } , 

{Jo X{x 1 {s)-q) J 

where ^(x, y) is the collection of piecewise C 1 paths x : [0, 1] — > V that satisfy 

(4.17) x(0) = (0, 0), x(l) = (a, y), and x'(s) G V for all s. 

The last condition ensures that 7 (x'(s)) is defined. The function A(-) in the 
definition of T q is the macroscopic rate defined by (2.12). Lastly, for q,x G R and 

t > 0, set 

(4.18) g«(x, t) = M{y : (x, y) G V, y) > t} . 

r 9 (a:, y) represents the macroscopic time it takes a £-type interface process to reach 
point (x, y). The point q marks the ^-coordinate of the defect column that (macro- 
scopically) has the slow rate A(0). The level curve of F q given by g q (- , t) represents 
the limiting interface of a ^-process, as stated in the next proposition. 
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Proposition 4.1. For all g,i6R and t > 0, 

(4.19) lim n~ l i\ nq \(nt) = g~ q (x,t) in probability. 

Remark 4-1- For the reader familiar with the proof of [Se3], let us point out that 
Prop. 4.1 is the analogue of Cor. 5.1 in [Se3], and T q corresponds to the limit in 
Prop. 5.1 in [Se3]. 

To prove Prop. 4.1 we follow Section 5 of [Se3] and switch to a last-passage 
representation. Lemma 4.2 below shows how the growth model of Section 2.1 
enters the picture, and justifies (2.7) as the definition of «(r). 

Define a lattice analogue of the interior of the wedge V by £ = G Z 2 : 

j > M > -j + 1}, with boundary DC = {(i,0) : % > 0} U {(i,-z) : i < 0}. For 
G £Ud£, let 

(4.20) L fc (i,j)=mf{t>0:^(t)>j} 

denote the time when £f reaches level j. The rules (4.9)-(4.11) give the boundary 
conditions 

(4.21) L k (i,j) = for (i,j)edC, 
and for G £ the equation 

(4.22) L k (i,j) = max{L fc (z - L k (i,j - l),L k (i + 1J- 1)} + ^ , 

where • is an exponential waiting time, independent of everything else. It repre- 
sents the time £ k waits to jump, after £ k and its neighbors £, k +1 have reached 
the positions that permit £ k to jump to j. By (4.10), I3 k _ k ^ has rate r, but for 
z 7^ — fe, /3, fc • has rate 1. 

The waiting time (3 k j cannot be read directly from T>i + k. One has to construct 
the evolution £ fc (-) up to the stopping time L = max{L fc (z — L k (i,j — 1), L k (i + 
l,j — 1)}, and then j3 k ■ is the waiting time to the next epoch in T>i + j~ after L. The 
last-passage representation entails switching probability spaces so that the waiting 
times f3 k j become the basic building blocks of the construction. We shall switch 
notation to keep the two constructions distinct. 

We now construct a last-passage growth model on £ that has a defect in the 
column {(m,j) : j G Z}, where m G Z is fixed. Let {r u : u G £} denote a collection 
of i.i.d. exponential rate 1 random variables. Define weights 



(4.23) w ™. 



1 , i ^ m 
1/r , i = m . 
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Given u = (u,v) G £, let II (u) denote the set of lattice paths n = {(0,1) = 
(«2,j2), • • • , (ip,j P ) = u} whose admissible steps satisfy 

(4.24) (i e ,j e ) - (ie-i,ji-i) = (1,0), (0,1), or (-1,1). 
Finally, define the passage times {T m (u)} by 

(4.25) T m (u) = for u G dC, 

and for u G £ as the maximal weighted sum of waiting times over admissible paths: 

(4.26) T m (u) = max V <Y V . 

7rGn(u) — ' 
v£7r 

Comparison of (4.21)-(4.22) and (4.25)-(4.26) shows that the passage time pro- 
cesses {L k (u) : u G £U<9£} and {T _fc (u) : u G £U<9£} have the same distribution. 
The superscript changes by a minus sign from L k (u) to T~ k (u) because L k (u) is 
the passage time for the process £ fc , but the slow column for this process is at —k, 
so we set m = —k in T m (u). 

The conclusion of this is that T~ k (i,j) has the same distribution as the time 
when reaches level j. So, as in [Se3], Prop. 4.1 will follow if we prove 

Proposition 4.2. For all q G R and (x,y) in the interior ofV, 

(4.27) lim n _1 T [n,7] ([nx], [ny]) = T q (x, y) in probability. 

Remark 4-2. Let us first explain the known homogeneous situation where r = 1 and 
A(0) = 1. Write ^ r=1 ) and T^ r=1 ) for the interface process defined by (4.9)-(4.11) 
and the passage time in (4.26) when r = 1 and there is no special column. In this 
case the limits in (4.19) and (4.27) are given by 

lim n~ 1 ^ r ~ 1 1 ' > (nt) = tq (x/t) = -( 1 — — ] for — t < x < t. and 

(4 . 28) r™ >*] ^ > ^ I > A\ t) ~ ~ 

lim n-^^dnx], [ny]) = 7o (x, y) = {^fx~Ty + ^/yf . 



The explicit values cannot be inferred directly from the path model, but indirectly 
via the connection with TASEP. Briefly, one computes the current fo(p) = p(l — p) 
from the known invariant distributions of TASEP. The coupling (4.8) implies that 
the limiting shape go for £ is the conjugate of /o, and one can derive go. The shape 
go is a level curve of the passage time 70, so one obtains the formula for 70 from 



(4.29) 



7o (a, 0o(aO) = 1, 



17 



and the homogeneity of 70. This type of argument is repeated in the examples in 
[Sel]. 

Proof of Proposition 4-2. Consider the possible maximizing macroscopic curves in 
(4.16). By the concavity of 70, of all the paths x(s) = (xi(s),X2(s)) in X(x,y), we 
only need to consider these two types: 

(4.30) x(s) = (sx, sy), < s < 1 

with value [except in the case x = q = which is covered by (4.32)] 

(4-31) / — f — 7T — ~^ds = n/o(x,y) , 



A(xi(s)-g) 



and 



^■(g,x 2 (si)), 0<s<si 
(4.32) x(s)={ (g,x 2 (si)) + ^;(0,x 2 (s 2 )-x 2 (si)), si < s < s 2 

(g,x 2 (s2)) + fEff (z -x 2 (s 2 )) , s 2 < s < 1, 



with value 



Jo 



70(X'( S) ) ^ 



(4.33) 7o A(xi(s) - q) 
= 7o(x(si)) + A(0)" 1 7o(x(s 2 ) - x( Sl )) + 7o ((x,y) - x(s 2 )) . 

In words: the path of type (4.30) is a single line segment from (0, 0) to (x, y). The 
path of type (4.32) first uses parameter interval [0, s\] to take a straight line path 
to the vertical line x = q, then spends interval [s\, s 2 ] on this line to take advantage 
of the slow rate A(0), and finally takes a straight line path to (x, y). 

Now consider the microscopic path problem (4.26). We need to establish the 
connection between it and the quantity A(0) = 4/«(r) defined in Section 2.1. 

Lemma 4.2. Set q = so that the special column goes through the origin. Then 
fory>0 

(4.34) lim n _1 T°(0, [ny]) = n{r)y in probability. 



Proof of Lemma The Lemma follows because the growth model of Sect. 2.1 
with a diagonal defect is the same as the one studied here when the columnar defect 
is at the origin. A simple mapping reveals this. First observe that the admissible 
step (0, 1) can be eliminated from (4.24), because each (0, l)-step in a path can be 
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replaced by a (1, 0)-step followed by a (—1, l)-step. This change adds a site to the 
path and hence increases its overall passage time. So we may assume that n(u) 
contains only paths that have admissible steps (1, 0) and (—1, 1). 

Consider the bijection ip : C — > N 2 given by = (i + For € N 2 , 

define random variables Yij = t^-i^j) and weights Wij = ^-i^jy Then Wij 
satisfies (2.4). For tt e 11(0, n), the image path i/j(tt) runs from (1,1) to (n,n), 
and has steps of two kind: (1,0) and (0, 1). Thus the map ip transforms T°(0, n) 
of (4.26) into T n n of (2.5). Now the lemma follows from the definition (2.7) of 
k(t). □ 

Return to the proof of Prop. 4.2. To first prove 
(4.35) lim inf n^T^Qnx], [ny]) > T q (x,y) , 

n^oo 

consider any macroscopic path x(-) of type (4.32). [We leave the easier type (4.30) 
to the reader.] Let 7r n be the microscopic path through the sites (0, 1), [nx(si)] = 
([nxi(si)], [nx2(s\)]), [nx(s2)], and ([nx], [ny]), constructed so that each of the three 
segments maximizes passage time between its endpoints. Then by (4.34) and the 
limit (4.28) for the homogeneous case, 



lim inf n-^^^jia;], [ny]) > lim inf n" 1 V J™ q W v 

n — >oo n — >oo < * 

VG7T„ 

> 7o(x(si)) + K(r)(x 2 (s 2 ) - x 2 (s 1 )) +io((x,y) - x(s 2 )) 
= the value of the path in (4.33). 

The last equality follows because A(0) _1 7o (x(s 2 ) — x(si)) = (l/4)«(r)7o(0, ^2(^2) — 
^2(^1)) = K>{r){x2{s2) — ^2(^1)). The reason there might not be equality in the 
last inequality in (4.36) is that an optimal path between, say, (0, 1) and [nx(si)] 
might actually take advantage of the [ng]-column and return a larger value than 
7o(x(si)). A similar argument for paths of type (4.30) justifies (4.35). 
Now for the complementary upper bound 

(4.37) limsupn-V^Qnz], [ny]) < T q (x,y) . 

Each macroscopic path x(-) in X(x, y) is contained in a fixed compact subset A of 
V. Choose 5 > so that 

(4.38) l7o(x) - 70 (y) I < £ for x, y G A such that |x - y| < 6, 
and then a partition 

= b < 61 < • • • < b s = y 
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of [0, y] with mesh max(6j +1 — bi) < 5. Let be the path of type (4.32) with 

x ( i> j \s 1 ) = b i <b j =xf j \s 2 ). 

Let us adopt the following generalization of the notation in (4.26): T m (u, v) 
denotes the maximal weighted sum over admissible paths from u to v. So T m (u) 
in (4.26) is the same as T m ((0, 1), u). 

Let n n be the maximizing microscopic path in (4.26) for u = ([nx], [ny]) and m = 
[nq]. The easy situation is when 7r n does not intersect the vertical column ([nq],j) 
that has the slow rate r. Then T^ nq \[nx], [ny]) equals the homogeneous passage 
time T^ r=l \[nx\, [ny]). If this happens infinitely often along the subsequence taken 
on the left-hand side of (4.37), then (4.37) follows from (4.28). 

Otherwise, pick indices k < I such that the path 7r n first touches the vertical 
column ([nq],j) in the range [nbk] < j < [nbk+i], and for the last time in the range 
[nb(\ < j < [nb£+i\. Then quite obviously 

T^([nx],[ny])= £ u^t v 

v6tt„ 

< T^([nq], [nb k+1 ])+T^(([nq], [nb k ]), ([nq], [nb e+1 ])) 
+ T^(([nq],[nb e ]),([nx],[ny])) 

x(t^= 1 )(H, [nb i+1 ])+T^(([nq], [nb,]), ([nq], [nb J+1 ])) 



+ T^(([nq],[nb J ]),([nx],[ny])) . 



< max 



Divide by n, let n — > oo, use the limits (4.28) and (4.34), and then (4.38) to get, 
withC = 1 + A(0) _1 , 



limsupn _1 T [n<zl ([n:r], [ny]) 



< max 70(9, b i+1 ) + K,(r)(b j+1 - h) + 7o(x - q,y - bj) 

< rn ax|7o(Q, bi) + A(0) _1 7o(0, bj - h) + -y (x - q,y - fy) J + Ce 

= maxjvalue of path x^' J ^(-)j + Ce 

i<j 

<T q (x,y) + Ce. 

This completes the proof of Proposition 4.2. □ 

Prop. 4.1 follows from Prop. 4.2 as Cor. 5.1 follows from Prop. 5.1 in [Se3]. 
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4-3. Hydrodynamic limit 

Using (4.12), (4.14), and (4.19) we can now prove that, in probability, 
(4.39) lim n^zVAnt) = v(x,t) = sup{v (q) - g~ q (x - q,t)\ . 

The argument is the one from eqn. (6.4) to (6.15) in [Se3], so we will not repeat it 
here. To complete the proof of (4.13) and thereby the proof of Thm. 2.2, we need 
to show that the limiting value v(x,t) defined above agrees with the desired limit 
v(x,t) defined by (2.14). 

4-4- Formula for v(x,t) 

By (4.16)-(4.18), the definition (4.39) of v(x,t) can be rewritten as 



(4.40) 



v(x, t) = sup < vo(q) — y : there exists a path x(-) e X{x — q, y) 

f 1 7o(x'(s)) 

such that / —, — — r- ds > t 

Jo \{xi(s)+q) 

Proposition 4.3. v(x,t) = v(x,t). 

Proof. The proof involves mapping the paths x(-) to the paths w(-) that appear in 
(2.14), and vice versa. 

Given a path x(-) G X(x — q, y) that appears in (4.40), define a new time variable 
T = t(s) by 

(4.41) r(.)= [' J° { * V) K ds. 



o X(x 1 (s) + q) 

Let the terminal r-time be t\ = r(l). From (4.40) we know that t\ > t. Let 
s = s(r) be the inverse time change. Define a path z : [0, t±] — > R by 

(4.42) z{r) = q + x 1 {s{r)). 
Then 

(4.43) z(0) = q , z{t\) = x , and z'(t) = x[(s(t))s'(t) . 

Differentiating (4.41), relation (4.29), and the homogeneity of 70 [means: 7o(cx, cy) = 
cio(x,y)} give 

(4.44) x' 2 (s) = r'(s)X(x 1 (s) + q)gJ A • 

\t'(s)X(x 1 {s) + q)J 
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Since y = #2(1), we can use this to compute 

• 1 f l 



V= f x' 2 (s)ds= f r\s)X(x 1 (s) +q)g ( ——^^——^ds 
Jo Jo \t'(s)X{x 1 {s) + q) J 

The only problem is that z(-) is defined on [0, t\] instead of on the possibly smaller 
interval [0,t]. Let w(a) = z(t\a/t) be a time change of z(-) defined for < a < t. 
Then, because sgo(x/s) is nondecreasing in s, change of variable r = {t\/t)a in the 
last integral above shows that 

v (q) ~y< v o (w(0)) - \{w{a))g Q {^^-^jda . 

Since x(-) was an arbitrary path inside the braces in (4.40), we have shown that 
v(x, t) < v(x, t). 

Conversely, take a path w : [0,£] — > R that appears in (2.14). Define a path 
x(cr) = (xi(a), x 2 (a)) for a G [0, 1] by 

xi(a) = w(at) — w(0) 

(4 - 45) ( \ r\( ( w fj^L^ 

Let q = w(Q) and y = 0:2(1). Then x(-) e X{x — q,y), provided x'(cr) e V [recall 
conditions (4.17)]. This follows because A and go are nonnegative functions, and 
because g(z) > —z for all z. Also, (4.29) and (4.45) give 

70 (xV)) _ 



A(xi(cr) + 5) 



so the integral condition inside the braces in (4.40) is satisfied. We conclude that 
x(-) is a path that appears in (4.40), and since 

v (w(0)) - [ \(w(s))g ( tt~7~tt ] ds = v (q)-y < v(x, t) , 
Jo \\(w(s))J 

we have v(x,t) < v(x,t). This completes the proof of Prop. 4.3. □ 
We have now proved Thm. 2.2. 



5. Proofs of the Corollaries 
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5.1. Proof of Corollary 2.1 

As a preliminary step for calculating macroscopic profiles from the variational 
formula (2.14), we optimize the integral term as a function of the initial point 
q = w(0). So let 



In terms of the limiting shapes of Prop. 4.1, I(x, t, q) = g~ q (x — q,t), so from the 
formulas below the reader can deduce explicit expressions for the limits in (4.19). 

By the convexity of go , it suffices to consider the following two types of paths in 
(5.1): either w(-) is a single linear segment from q = w(0) to w(t) = x; or it consists 
of a linear segment from q = w(0) to w(si) = 0, a constant segment w(s) = for 
si < s < S2, and a linear segment from w(s2) = to w(t) = x. Only calculus is 
involved in finding the optimal paths, so we skip the details and present a summary 
of the results. Abbreviate 



(5.1) 




B=y/1- A(0) . 



Five different ranges of the variables x and q appear. 



(5.2a) 



x\ > Bt ; 



< x < Bt , and 



(5.2b) 




-Bt < x < , and 



(5.2c) 




< x < Bt , and 



(5.2d) 



2 

x - Bt < q < (VBt - y/x) ; 



—Bt < x < , and 



(5.2e) 
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These are the optimal values: 

Cases (5.2a-c): I(x,t,q) = tgo((x — q)/t) and the optimal path is w(s) = q + 
(s/t)(x-q). 

Cases (5.2d-e): 

I(x,t,q) Mg (-Bq/\q\) + - M + M) \(0)g (0) + ^<7o (Bx/\x\ ) 

=f (M + kl)-^=^ + |A(o) 

and the optimal path is 

{q — sq/si , < s < S\ 

, Si < S < S2 

(s - s 2 )x/(l - s 2 ) , s 2 <s<l, 

with 

Si = \q\/B and s 2 = t — \x\/B . 

Proof of Corollary 2.1 is now reduced to finding 

(5.3) v(x,t) = sup{v (q) - I(x,t,q)} 

with vo(q) = pq. We skip the calculus details. 

5.2. Proof of Corollary 2.2 

Assume now that p* < 1/2, in other words, that the slow bond disturbs the 
hydrodynamic profiles. First check from (5.3) that all the profiles admitted by the 
restrictions stated in Cor. 2.2 are in fact invariant. This contains the following cases: 
constants at p* and at 1 — p*; a piecewise constant profile with a single entropy 
shock [jump from p* to 1 — p*] anywhere in R; a piecewise constant profile with a 
single non-entropy shock [jump from 1 — p* to p*] at x = 0; and a piecewise constant 
profile with a non-entropy shock at x = 0, and an entropy shock in (— oo, 0), or in 
(0, oo), or in both. 

To prove that these are the only possible invariant profiles, let po be a profile 
such that p* < po(x) < 1 — p* and po is invariant under (5.3). Define vq by (2.11). 
The invariance means that v x (x,t) = po(x) a.e., and so 

(5.4) v(x,t)= [ p o (y)dy + v(0,t) = v o (x) + v(0,t). 

Jo 
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Formula (5.3) operates like a semigroup, and this with (5.4) implies that v(0, t) = Ct 
for some constant C. To determine C, we do a comparison. Let 

vl(x) = p*xl{x < 0} + (1 - p*)xl{x > 0} , v e (x, t) = v%(x) - £A(0)/4 

and 

v£(x) = (1 - p*)xl{x < 0} + p*xl{x > 0} , v n (x, t) = v%(x) - £A(0)/4 

denote the evolution of the entropy and the non-entropy shock at the origin. The 
bounds v n (x,t) < v(x,t) < v e (x,t) are valid because they are valid at time t = 
and preserved by (5.3). Consequently 

(5.5) v(x,t) = v (x) -U(0)/4. 

Lemma 5.1. There cannot exist an x ^ 0, and e, 5 > such that this holds: 

v (x) > v (q) + (p* + S)(x- q) for qE [x-e,x] 
and vq(x) > vo(q) — (1 — p* — 8)(q — x) for q e [x, x + e] . 



Proof. Suppose such x, s, 5 exist. Pick t > small enough so that t < e and 
Bt <\x\. Then, by (5.2), (5.3) becomes 

(5.7) v(x,t) = sup{v (q) - tg ((x - q)/t)} . 

gen 

We shall show that formula (5.7) gives something strictly smaller than (5.5), and 
this contradiction makes (5.6) impossible. 

In (5.7) it suffices to consider q G [x — £, x + 1], by observing from (2.13) that go 
has constant slopes to the left of —1 and to the right of 1. For q G [x, x + 1] write 

vo(q) - t9o((x - q)/t) < v Q (x) + (1 - p* - S)(q - x) - tg ((x - q)/t) 

= t*(z)-f{(l-p*- *)£ + <*,(£)}, 

where ^ = (x — q)/t may vary freely in [—1,0]. By the duality of g and the TASEP 
current fo(p) = p(l — p), the expression in braces is bounded below by fo(l — p* —6) 
= fo(p* + 5). This and a similar argument for q G [x — t, x] give 

v(x,t)<v (x)-tf (p* + S). 

Since A(0)/4 = fo(p*), the above bound is strictly less than (5.5), provided 6 is 
chosen small enough to have p* + 6 < 1/2. □ 
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Now we can prove that p takes only the values {p*, 1 — p*}, up to Lebesgue 
null sets. For suppose the set A = {x : p* + 5o < po(x) < 1 — p* — 5o} has positive 
Lebesgue measure for some 5q > 0. Then A has a density point x. (For a definition, 
see for example p. 107 in [WZ].) We show that (5.6) holds at x. Let q < x. Since 

Po(x) > p*, 

v (x)-v (q)= / p (y)dy+ / Po(y)dy 

J[q,x]CiA J[q,x]nA c 

> p*(x - q) + 5 ■ m([q, i]nA), 

where we wrote m for Lebesgue measure. Since a; is a density point, m([g, x] C\A) > 
(x — q)/2 if q is close enough to x. This checks the first part of (5.6), and the other 
part is similar. 

Assuming that po is piecewise continuous, we now conclude that it is piecewise 
constant with values p* and 1 — p*. To prove Cor. 2.2, it remains to observe that 
Lemma 5.1 prevents a jump from 1 — p* to p* everywhere else except at x = 0. 



5. 3. Proof of Corollary 2. 3 



Let the initial distribution of the process be the i.i.d. product measure cto with 
density ao{i]i = 1} = p, with p outside the disturbed range (p*, 1 — p*). Let a t be 
the distribution of the process at time t. Let p be a limit point of the time averages 
of a t , so for some sequence tk S oo, 



1 '* 



p = lim — / a s ds 

k^oo tk Jo 

in the weak sense on the compact state space {0, 1} Z . Such a limit point exists by 
compactness, and is then automatically invariant for the process. 

Write E 010 for expectation under the path measure of the process started with 
distribution ao- For i ^ we have the equation 

(5.8) E ao [Ji(t)] = [ a s { Vi = 1, Vl+1 = 0}ds 

Jo 

because Ji(t) [= the number of jumps from site i up to time t] increases by 1 at 
rate 1 when the event inside the braces holds. Thus 

1 f tk 

p{r]i = 1, rji+i = 0} = lim — / a s {r)i = 1, rj i+1 = 0}ds 

fc^oo tk Jo 

(5 ' 9) = lim ^E^lUtk)] 

= p(l-p). 
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The last equality is from (2.17) and Case 1 of Cor. 2.1, with initial function Vq(x) = 
px. Ji(t) is bounded by a Poisson(t) random variable, hence there is uniform 
integrability to justify the limit above. If i = 0, we have to multiply the right-hand 
side of (5.8) by the factor r, so the limit on the last line of (5.9) is multiplied by 
r . This proves (2.21) for p. 

Suppose the initial distribution «o has a macroscopic profile po in the sense of 
(2.16) such that p* < po(x) < 1 - p*. Then the limit in (5.9) is valid with p = p*. 
In particular, this is the case for the cvq with non-entropy shock described in the 
second paragraph of Cor. 2.3. 
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